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^-H ■ Abstract 
> ■ 

, The neutrino spin-flavor conversion of Ve <-> which is induced by the inter- 

' action of the neutrino magnetic moment and magnetic fields in the collapse-driven 

. supernova is investigated in detail. We calculate the conversion probability by us- 

\^ I ing the latest precollapse models of Woosley and Weaver (1995), and also those of 

0^ ' Nomono and Hashimoto (1988), changing the stellar mass and metallicity in or- 

, der to estimate the effect of the astrophysical uncertainties. Contour maps of the 

conversion probability are given for all the models as a function of neutrino mass 
^ ' squared difference over neutrino energy (Am^/iJ^) and the neutrino magnetic mo- 

^ , ment times magnetic fields {/ij^B). The expected deformation of the Ue spectrum can 

^ ' be clearly seen from these maps, and some qualitative features which will be useful 

^ ■ in the future observation are summarized. It is shown that in the solar metallicity 

, models some observational effects are expected with Am^ = 10~^-10~^ [eV^] and 

' /i^ > 10~^^(10^G/-Bo) [ij-b], where Bo is the strength of the magnetic fields at the 

^ . surface of the iron core, and Bo = 10^ [Gauss] is a reasonable value which is inferred 

from the observed magnetic fields in white dwarfs. We also find that although the 
dependence on the stellar models or stellar mass is not so large, the metallicity of 
precollapse stars has considerable effects on this conversion. In lower metallicity 
stars, strong precession between and occurs with smaU Am'^/E„ (< 10"** 
[eV^MeV]), and conversion probability changes periodically with Bq. Such effects 
may be seen in a supernova in the Large or Small Magellanic Clouds, and should be 
taken into account when one considers an upper bound on fj,i, from the SN1987A 
data. 
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I. INTRODUCTION 



Copious neutrino emission from collapse-driven supernovae attracts significant attention because 
it provides rich information not only on the mechanism of supernovae but also on the neutrino 
physics through a number of events captured in some underground neutrino detectors, such as 
the Super-Kamiokande (SK) j^. The most noteworthy subject on the nature of neutrinos is the 
mass of neutrinos and oscillations between different flavors induced by the mass difference. How- 
ever, the ordinary matter oscillation (well known as the MSW effect |Q) has its effect only on 
neutrinos but not on antineutrinos under the direct mass hierarchy, and the vacuum oscillation 
is not observable unless the mixing angle is unnaturally large compared with that of the quark 
sector. In this case, electron antineutrinos (S'e's), which is the most detectable in a water Cerenkov 
detector, do not undergo any oscillation. One of some possibilities that S-e's would oscillate or be 
converted into other species of neutrinos is the neutrino magnetic moment. If the neutrinos have 
a nonvanishing magnetic moment, it couples the left- and right-handed neutrinos, and interaction 
with sufficiently strong magnetic fields induces the precession between neutrinos with different 
chiralities in the inner region of the collapse-driven supernova |^,^ . In general, non-diagonal ele- 
ments of the magnetic moment matrix are possible, and neutrinos can be changed into different 
flavors by this flavor changing moment ||^. Furthermore, with the additional effect of the coherent 
forward scattering by dense matter in the collapsing star, neutrinos can be resonantly converted]^ 
into neutrinos with different chiralities by the mechanism similar to the MSW effect. This 

resonant spin-flavor conversion induced by the neutrino magnetic moment may drastically deform 
the spectrum of electron antineutrinos (Pe's) in the water Cerenkov detectors. The earlier publica- 
tions have shown that in the future experiment this effect will be observable with inner magnetic 
fields of some reasonable strength, if there is a magnetic moment a little smaller than the current 
astrophysical upper limits from the argument of the stellar cooling due to the plasmon decayQ: 
< 10~^^-10^^'^fj,B, where /is is the Bohr magneton Jl2| , p^ . The magnetic moment of neu- 
trinos in the standard electro-weak theory with small neutrino masses is very small due to the 
chirality suppression; for example, the standard SU(2)lxU(1) model with a singlet right-handed 
neutrino gives /ii, ^ 3 x 10~^^(mi,/leV)^B, far below the experimental/astrophysical upper bounds 
4|-16|. However, some particle-physics models [^,^ have been proposed in order to give a 



large magnetic moment of ^ lO^^^fis which would explain p9|"|21|| the anticorrelation between 
the time variability of the solar neutrino flux and the sun spot numbers suggested in the '^''Cl 
experiment ||2^ . (The anticorrelation in the CI experiment, however, has not yet been statistically 
settled.) Therefore, the influence of a large magnetic moment on various physical or astrophysical 
phenomena including collapse-driven supernovae deserves more detailed investigation. Discovery 
of a large magnetic moment of neutrinos indicates that there exist interactions which violate the 
chirality conservation beyond the standard theory. 

In this paper, the resonant spin-flavor conversion between right-handed S'e's and left-handed 



^The precession itself is suppressed by the matter potential. 

^This constraint refers to the norm of the neutrino magnetic moment matrix, (X/i j Im^j P)^''^! i-^-i ^^is 
includes the flavor changing moment. 
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mu or tau neutrinos (i'^ or z^r's) is studied assuming that the neutrino is the Majorana particle. 
In general, the matter potential suppresses the interaction of the magnetic moment and magnetic 
fields because of the generated difference of the energy levels. However, Athar et al. pointed out 
that the resonant conversion of this mode {ve ^ v^i) occurs quite effectively in the region above the 
iron core and below the hydrogen envelope of collapsing stars, namely, in the 0+Si, 0+Ne+Mg, 
0+C, and He layers (hereafter referred to 'the isotopically neutral region'). The reason is that 
the effective matter potential for the De i^^i mode is given in the form proportional to the value 
of (1 — 2ye), where Yg is the electron number fraction per nucleon, and is very close to 0.5 in 
this region (typically, (1 — 2Ye) ^ 10~'*-10~^); the matter effect is therefore strongly suppressed 
compared with the magnetic interaction, and the adiabaticity condition becomes considerably less 
stringent. Athar et al. have shown that assuming ^i,{= /ip^i/^) ^ lO^^^/is, this resonant 
conversion would occur with some reasonable assumptions about magnetic fields in a star. We also 
consider around this valuc^. 

In order to judge the deformation of an observed spectrum as the evidence of the existence 
of the neutrino magnetic moment, it is necessary that the conversion probability is calculated 
with high accuracy in a wide range of some parameters such as the mass of neutrinos or the 
magnetic fields. However, only rough estimates or demonstrations in some cases are given in the 
earlier publications and the relation between the shape of the deformed spectrum and the physical 
parameters has not yet been clarified. Therefore we make the contour maps of the conversion 
probability of D,, <-> for some used models of precollapse stars as a function of the two parameters: 
t^rr? jEy and n^Bo, where Am^ is the neutrino mass squared difference, E^, the neutrino energy, 
and Bq the magnetic field at the surface of the iron core. The expected observational effects can 
be clearly understood by these maps. Some examples of spectral deformation are also calculated 
and qualitative features which will be useful for the future observation are summarized by using 
these maps. 

It is apparent that the deviation of the value of Ye from 0.5 in the isotopically neutral region is 
quite important, and this value is strongly dependent on the isotopic composition. Since almost 
all nuclei in the isotopically neutral region are symmetric in the number of neutrons and protons, 
this deviation is determined by rarely existent nuclei and the accurate estimate of this deviation is 
quite difficult. Therefore, the astrophysical uncertainty in (1 — 2Ye) should be discussed. We use 
the latest 15 and 25 Mq precollapse models of Woosley and Weaver (hereafter WW) which 
include no less than 200 isotopes. Such a large number of isotopes have never been used previously 
in the calculation of (1 — 2Ye). It is also expected that this value strongly depends on the stellar 
metallicity, and hence we use the WW models with the two different metallicities: the solar and 
zero metallicity, and the metallicity effect is investigated. Also the models of 4 and 8 Mq helium 
core of Nomoto & Hashimoto (hereafter NH) are used, which correspond approximately to 15 
and 25 Mq main sequence stars, and the model dependence of (1 — 2Ye) is discussed. 



^There is a further stringent constraint on the transition magnetic moment of massive neutrinos from 
observation of the 21-cm (hyperfine) radiation from neutral hydrogen gas in external galaxies: fi^, < 
1.7 X 10~^^ for the neutrino masses above 30 eV [ p3|j2^ . However, we consider only the range of Am^ < 1 
[eV'^] and this upper bound does not constrain our analysis. 
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We consider Am^ smaller than about 1 [eV^], therefore the resonance occurs above the surface of 
the iron core. The resonance in the iron core and its implications on the dynamics of the supernova 
considered in a recent preprint p7| are not discussed here. The global structure of magnetic field is 
assumed to be a dipole moment, and the strength of the magnetic field is normalized at the surface 
of the iron core with the values of 10®-10^° [Gauss], which are inferred from the observation of the 
magnetic fields on the surface of white dwarfs. Throughout this paper, we consider the conversion 
between two generations for simplicity. Because and i^r can be regarded as identical particles 
in the collapse-driven supernova, our results also apply to the conversion of De v-r- Derivation of 
the equation which describes the propagation of neutrinos and evolution of conversion probability 
is given in section pi and the profile of the effective matter potential and magnetic fields are 
given in section |ll| by using the precollapse models of massive stars. Qualitative features of 
the conversion are also discussed in this section. Numerical results are given in section 1^, and 
spectral deformation is also discussed. Discussion and conclusions are given in sections ^ and [v| , 
respectively. 



II. FORMULATIONS 

The interaction of the magnetic moment of neutrinos and magnetic fields is described as 

< {l^z)R\Hint\{t^j)L f^ijB± , (1) 

where fiij is the magnetic moment matrix, B± the magnetic field transverse to the direction of 
propagation, {i')r and the right- and left-handed neutrinos, respectively, and i and j denote 
the flavor eigenstate of neutrinos, i.e., e, fi, and r. The magnetic moment interacts only with 
transverse magnetic fields. If neutrinos are the Dirac particles, right-handed neutrinos and left- 
handed antineutrinos do not interact with matter and therefore undetectable. The conversion into 
these sterile neutrinos due to the magnetic moment suffers strong constraints from the observation 
of neutrinos from SN1987A by the Kamiokande II |^ and 1MB |2^, and also from the argument 
on energy transportation in the collapse-driven supernova ||^, p^ , ^0| -^ . On the other hand, if 
neutrinos are the Majorana particles, as assumed in this paper, vrs are antineutrinos and interact 
with matter, and the constraint becomes considerably weak. The diagonal magnetic moment is 
forbidden for the Majorana neutrinos, and therefore only the conversion between different flavors 
is possible, e.g., (i^e)i? {i^ii,t)l- As mentioned in introduction, we investigate this mode because 
the conversion of this mode occurs quite effectively in the isotopically neutral region and also u^s 
are most easily detected in the water Cerenkov detectors. 

In dense matter of the collapsing stars, the coherent forward scattering by matter leads to 
the effective potential for neutrinos, and this potential for each type of neutrinos is determined 
according to the weak interaction theory. The potential due to scattering with electrons is given 
as (including both the charged- and neutral-current interactions) 

V• = ±^/2 Gf (±^ + 2sin2 6'vF) , (2) 

where Ue is the number density of electrons, Gp the Fermi coupling constant, and 9w the Weinberg 
angle. The ± sign in the parentheses refers to {+) and Vf^^r (— ), and that in front to v {+) 
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and />(—). In the ordinary flavor oscillation (v^ ^ ^ii.t), the effective potential is only due to 
the charged-current scattering by electrons because the effect of neutral-current interactions is the 
same for all flavors. However, we have to consider the neutral-current interaction in the conversion 
of v and z/, because of the opposite signs of the potential. Therefore the neutral-current scattering 
by nucleons should also be included, that is 

y = ±V2 Gf -2sin2 6'H/) ripT V2 Gf ^ n„ , (3) 

where Up is the number density of protons, n„ that of neutrons. The ± or =p signs refer to u 
(upper) and v (lower) for all three flavors of neutrinos. We do not have to consider the form factor 
of nuclei because the relevant interaction is forward scattering and there is no momentum transfer. 
The isotopically neutral region is far beyond the neutrino sphere and neutrinos go out freely in this 
region; hence we do not have to consider the neutrino-neutrino scattering. By using the charge 
neutrality, the difference of the potentials for i/g's and f^'s (or z^r's) which we are interested in is 
as follows: 

Ay = Vu^ - = V2GFp/mN{l - 2Ye) , (4) 

where p is the density, niN the mass of nucleons, and ~ np/{np + n„). Now the time evolution 
of the mixed state of and is described by the following Schrodinger equation: 




p^B^ 
^^B_L Ai? 




(5) 



and /S.H is defined as: 



Am^ 

AF= -— cos2e'- AF , (6) 

where -B^ is the energy of neutrinos, Am^ = ra\ — m\, 9 the angle of the vacuum generation mixing, 
and r the radius from the center of the star. Here we consider only and u^, but this equation 
is actually a truncation of the original 4-componcnt (i^g, v^^ and antincutrinos) equation (see ref. 
1^). The neutrino masses, mi and m2 are those in the mass eigenstates (m2 > mi). The direct 
mass hierarchy is assumed here and therefore Am^ is positive. The other terms have their standard 
meanings and the units of c = ?i = 1 are used. Also note that we can subtract an arbitrary constant 
times the unit matrix from the Hamiltonian, which does not affect the probability amplitudes. In 
the MSW flavor oscillation, there appears the term of generation mixing, Am^ sin20/4i?y, in the 
off-diagonal elements of the Hamiltonian; however, this term does not appear in this spin-flavor 
conversion between neutrinos and antincutrinos. In the following, pi, and cos 20 are set to be 
l{)^^^p,B and 1, respectively, and the scaling of B or Arn^ with respect to other values of p^ or 
cos 29 is obvious. The resonant spin-flavor conversion occurs when the difference of the diagonal 
elements in the Hamiltonian vanishes, and hence the resonance condition is given as AiJ = 0. By 
using this equation, the probability of conversion can be calculated provided that p{r), Ye{r)^ and 
B±{r) are known. 
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III. ASTROPHYSICAL ASPECTS 

A. Effective Matter Potential 

In this section, we consider the effective matter potential in the isotopically neutral region. The 
value of (1 — 2Ye) which we are interested in is easily calculated as: 




where Zi, Ai, and Xi are the atomic number, mass number, and the mass fraction of the i-th 
isotope, respectively, and the subscript i runs over all isotopes with 2Z ^ A. In order to get 
this value and the density profiles, the precollapse models of massive stars of Woosley &: Weaver 
(WW) and Nomoto & Hashimoto (NH) are used. We assume that the dynamical effect 
can be ignored within the time scale of the neutrino emission, and hence use the above static 
models. The mass and radius of the helium core of a 15Mq main sequence star is ~ 4Mq and 
~ li?0, respectively, and its free-fall time scale, {^/Gp)~^ is ~ 10^-10'^ [sec], which is longer than 
the neutrino emission time scale (at most a few tens of seconds). It takes about several tens of 
seconds for the shock wave generated at the core bounce to reach the hydrogen envelope , and 
the inner region of the isotopically neutral region may be disturbed by the shock wave. We will 
discuss about this in section |^. The calculation of the WW models of 15 and 25 Mq (hereafter 
WW15 and WW25, respectively) includes 200 isotopes, up to ''^Ge. Although the network of 19 
isotopes is used for energy generation up to the end of oxygen burning, the network of 200 isotopes 
is updated in each cycle and mixed using the same diffusion coefficients. The NH models are 4 
and 8 Mq helium cores (hereafter NH4 and NHS, respectively) corresponding approximately to 15 
and 25 Mq main sequence stars. Their calculation includes 30 isotopes up to the end of oxygen 
burning, which are also used for the energy generation. The WW and NH models use the different 
reaction rates of "'^^C(a, 7)^^0, and the treatment of convection is also different. As for the WW 
models, we use the models with two different metallicities; the solar and zero metallicity, and the 
metal abundance of the NH models is that of the Sun. 

By using the data of composition as well as the density profile of the solar metallicity WW 
models (WW15S and WW25S, where 'S' denotes the solar metaUicity), \AV\ in the WW15S and 
WW25S models are depicted in Figs. ^ and ^ respectively, by the thick solid lines as a function of 
the radius from the center of the star. Also shown by the dashed lines is \AH\ when /S.m? /Ei, is 
10-* and 10~^ [eV^/MeV]. The dominantly existent nucle i are also indicated for each layer in the 
top of these figures. In the neutronized iron core, Ye is smaller than 0.5, and AV^ is positive and 
much larger than /S.m? /Ei, unless /S.m? /Ei, is larger than 10^^ [eV^/MeV]. We consider the range of 
lS.m? I El, below this value, and hence the resonance does not occur in the iron core. Above the iron 
core, i.e., in the isotopically neutral region, Ye becomes quite close to 0.5 (still Ye < 0.5) and AV 
is strongly suppressed, typically by a factor of ^ 10^"^ in solar metallicity stars, and the term /x^S 
becomes more effective. This suppression continues to the end of the isotopically neutral region, 
namely, just below the hydrogen envelope. The isotopically neutral region is roughly divided into 
the four layers: O+Si, O+Ne+Mg, O+C, and He layer, from inner to outer region. The values of 
(1 — 2Ye) and some nuclei which are relevant to the deviation of Ye from 0.5 are tabulated in Table 
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H for each layer and for the six precoUapse models used in this paper. For the solar metallicity 
models, (1 — 2Ye) is determined mainly by the isotopes such as ^^Ne, "^^'"^^Mg, ^^Al, ■^''S, '^'^Ar, and 
so on. 

The resonance occurs when AV becomes smaller than Arri^ /2E^, and after the resonance (above 
the resonance layer) AH becomes constant with radius because AV is negligibly small. If the 
strength of the magnetic field is sufficiently strong for the satisfaction of the adiabaticity condition 
at the resonance layer, the neutrinos are resonantly converted into the other helicity state. The 
magnetic fields and the adiabaticity condition are discussed in the following subsection. The 
resonance layer is in the isotopically neutral region if Am'^/E^ is in the range of roughly 10^^°- 
10^^ [eV^/MeV] (slightly dependent on the stellar models), and the resonance layer moves inward 
with increasing Am? / E^. If Am? /En is smaller than 10"^" [eV^/MeV], the mass term has no effect 
on AH in this region of the solar metallicity models, and the resonance occurs at the boundary 
between the helium layer and the hydrogen envelope due to the change of the sign of (1 — 2Ye). In 
contrast with the flavor oscillation, the matter potential changes its sign by itself and the resonance 
can occur without the mass term. Am? jEy. However, as explained in the next subsection, if a 
dipole moment is assumed as the global structure of the magnetic fields, it seems difficult that B is 
strong enough to satisfy the adiabaticity condition at this boundary in the solar metallicity models. 
In the hydrogen envelope. Ye is about 0.8 and the suppression of (1 — 2Y^) does not work any more. 
We can see that most of the qualitative features are the same for the two models: WW15S and 
WW25S, and the dependence on the stellar masses is rather small. Note that our result gives 1-2 
orders of magnitude larger than that in the earlier calculation by Athar et al. in which 
the older 15 Mq Woosley & Weaver model fsj,^ is used. It is probably because our calculation 
of (1 — 2Ye) includes the larger network of isotopes used in the latest WW models. 

In Figs. ^ and ^ we show the same with Figs. Q and ||, but for the Nomoto & Hashimoto 
models. It can be seen that the profiles of AV of the WW and NH models are not so different, 
and the model dependence seems rather small. However, the situation is drastically changed when 
we consider the effect of different metallicities. Figs. || and ^ are the same with Figs. |l| and ||, but 
for the zero metallicity WW models: WW15Z and WW25Z. ('Z' denotes the zero metallicity.) In 
the 0+Si and 0+Ne+Mg layers, (1 — 2ye) is smaller than that of the solar metallicity models by 
about 1 order of magnitude, and in the 0+C and He layers, (1 — 2ye) is further strongly suppressed 
(4-6 orders of magnitudes) because of the lack of the heavy nuclei which cause the deviation of 
Ye from 0.5 (Table ||). In consequence, the metallicity effect becomes especially important when 
Aw? I El, is smaller than ~ 10~^ [eV^/MeV]. How this effect changes the profile of the conversion 
probability will be discussed in more detail in section 



B. Magnetic Fields 

Let us consider the magnetic fields in the isotopically neutral region. In the earlier publication 
, the strength of magnetic fields was normalized at the surface of the newly born neutron star 
(r ^ 10 km), but it is unlikely that the magnetic fields of a nascent neutron star have some effects 
on the far outer region, such as the isotopically neutral region, within the short time scale of 
the neutrino burst. The magnetic fields should be normalized by the fields which are static and 
existent before the core collapse. The strength of such magnetic fields above the surface of the iron 
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core may be inferred from that observed on the surface of white dwarfs, because the iron core of 
giant stars is similar to white dwarfs in the point that both are sustained against the gravitational 
collapse by the degenerate pressure of electrons. The observations of the magnetic fields in white 
dwarfs show that the strength spreads in a wide range of 10^-10^ Gauss [|6). Taking account of 
the possibility of the decay of magnetic fields in white dwarfs, it is not unnatural to consider the 
magnetic fields up to 10^" Gauss at the surface of the iron core. As for the global structure of the 
fields, although the optimistic estimate oi B oc is sometimes discussed from the argument of the 
flux freezing, a magnetic dipole is natural as static and global fields; we hence assume such fields 
in this paper. Therefore the off-diagonal element of the Hamiltonian in Eq.(||), iiuB± becomes 
fi^Bo{rQ/r)^ siiiQ, where Bo is the strength of the magnetic field at the equator on the iron core 
surface, rg the radius of the iron core, and Q the angle between the pole of the magnetic dipole 
and the direction of neutrino propagation. If the magnetic field is normalized at the surface of the 
neutron star, the radial dependence of a dipole (ex r~^) gives too small field in the isotopically 
neutral region, but the normalization at the surface of the iron core inferred from the observations 
of white dwarfs makes it possible that the magnetic field is sufficiently strong in the isotopically 
neutral region under the condition of a global dipole moment. The lines of fj,,^B [eV] are shown in 
Figs. |^-||for Bq = 10* and 10^" [Gauss], assuming fi^, = 10~^^ fiB- The strength of magnetic fields 
is also discussed from the argument of energetics. The energy density of the maximum strength of 
magnetic fields should at most be the same order of magnitudes with that of the thermal plasma 
in the star. Let us define the magnetic fields Bth, whose energy density is the same with that of 
the gas in the star: 

LBf,^^-^kT , (8) 
Stt 2 limp ' ^ ' 

where p is the density, T temperature, k the Bolzman constant, and fl the mean molecular weight. 
The line of iiyBth is depicted in Figures assuming fi. — 1 and = lO^^^/is, and it can be 
seen that the magnetic fields up to Bq ~ 10^" [Gauss] are far below Bth and therefore natural from 
the view point of energetics. 

If there is no matter potential, the complete precession of vb. ^ vl occurs; however, the pre- 
cession amplitude is suppressed by the matter potential. The precession amplitude is given in the 
form 0Jl|]: 

In the neutronized iron core, AiJ is much larger than /i^B even when Bq ^ VS^^ [Gauss], and the 
precession below the surface of the iron core can be completely neglected. (In other words, we 
can start the calculation with the pure neutrino states from the iron core surface, with i?o ^ 10^'^ 
[Gauss].) Above the iron core, i.e., in the isotopically neutral region in the solar metallicity models, 
[LyB is still much lower than ISV (or A_ff), except at the resonance layer or the boundary of the 
helium layer and the hydrogen envelope, as shown in Figs. Therefore, the precession does not 
occur in the solar metallicity stars. However, if the strength of the magnetic fields is strong enough 
to satisfy the adiabaticity condition at the resonance layer, neutrinos are resonantly converted into 
other types of neutrinos. The adiabaticity condition is satisfied when the precession length at the 
resonance layer, {^i^B)^^, is shorter than the thickness of the resonance layer, i.e.. 
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1 1 

2 

(at the resonance), (10) 
since the thickness of the resonance layer, Ari.os7 is given as 

Ar^os = fJ.„B 

Note that the suppression of (1 — 2Ye) in the isotopically neutral region makes the adiabaticity 
condition well satisfied because it reduces the right hand side of Eq.(|lO|) by a factor of (1 — 2ye)^^^- 
In order to show how this condition is satisfied, \d{AV)/dr\^^^ is shown by the thin solid lines 
in Figs. If /ii/B is (roughly) larger than \d{AV)/dr\^^'^ at the resonance layer {AH = 0), 

the adiabaticity condition is satisfied and De's and f^'s are mutually converted. In both the WW 
and NH models with the solar metallicity, the region where this condition is satisfied appears with 
Bq ^ 10^" [Gauss]. Because the slope of \d{AV)/dr\^/^ is fiatter than that of /i^S, this condition 
is satisfied better in the inner region of the star, in other words, with large values of Am? / E^, 
in the solar metallicity models. When Aw? jEy is smaller than ^ 10~^° [eV^/MeV] and the 
resonance layer lies at the boundary of the helium layer and the hydrogen envelope, unnaturally 
strong magnetic fields are necessary for satisfaction of the adiabaticity condition. However, in the 
zero metallicity stars, because the value of (1 — 2Yf,) is very strongly suppressed in the 0+C and 
He layers, [lyB becomes much larger than AH and hence the strong precession between different 
chiralities occurs with small Am^/Ey (Figs. || and||). Since the adiabaticity may be broken at 
the quite large jump of the matter potential at the boundary of the helium layer and the hydrogen 
envelope, the detailed calculation is necessary for the conversion probability when Am? /E^ < 10~^ 
[eV^/MeV] . It is apparent that the conversion probability in zero metallicity stars will be completely 
different from the solar metallicity stars. Now all of the qualitative features of the conversion can 
be understood from Figs. and the results of final conversion probability obtained by solving 
the evolution equation (Eq. H) numerically are given in the following section. 



B> 



d{AH) 
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d{AV) 


dr 




dr 



d{AH) 
dr 



(11) 



IV. RESULTS 

A. Conversion Probability Maps 

In this section the contour maps of the conversion probability {D^ <-> v^j) are given for all the 
models used in this paper as a function of Aw? j E^, and B at the surface of the neutronized iron core 
(Bq)- Before we proceed to contour maps, the evolution of conversion probability in the isotopically 
neutral region along the trajectory of neutrinos is shown for some cases as a demonstration. Fig. |^ 
shows the conversion probability as a function of radius from the center of the star using the model 
NH4 for some values of Bo, with Am^jE^ = 10-4[eVVMeV], y.^, = lO^^Vs, and cos 26* = 1. The 
resonance layer lies at r '--^ 5 x lO~^i?0 in the 0+Si layer and its location is never changed by 
strength of the magnetic fields. We can see in this figure that the conversion probability becomes 
larger with increasing strength of magnetic fields, and the complete conversion occurs with the 
magnetic fields strong enough (Bq 5 x 10^ [Gauss], in this case) to satisfy the adiabaticity 
condition (see also Figure ||). Fig. || is the same with Fig. |^, but Aw? jE^ is 10^^[eV^/MeV] and 
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the resonance layer is hence in more outer region at r ~ 1.5 x lO^^i?© (0+Ne+Mg layer). As 
shown in these figures, the necessary Bq for the complete conversion becomes larger with decreasing 
Am^ /El, in the solar metallicity models, because the adiabaticity condition is well satisfied with 
larger Am^/E^, as discussed in the previous section. In both the figures, the conversion probability 
jumps up a little at the radius of about 2.5 x lO~^i?0, because this radius corresponds to the surface 
of the iron core and the value of (1 — 2Yi,) drops quite suddenly here. 

Now we calculate the contour maps of the conversion probability for the solar metallicity models, 
in the region of Arn^jEy = 10-S-10-i[eVVMeV] and Bq = 10«-10i° [Gauss], and the results 
are given in Figs. ||-|l| (for WW15S, WW25S, NH4, and NHS, respectively). In the region of 
Avr? jEy < 10~®[eV^/MeV] or Bq < 10® [Gauss], the conversion does not occur because of too 
weak magnetic fields. Magnetic fields stronger than 10^" [Gauss] induce the precession below the 
surface of the iron core which cannot be ignored, and Am? /Eiy larger than 10~^[eV^/MeV] leads 
to the resonance below the surface of the iron core. In this paper, we consider the parameter region 
in which the conversion or precession below the iron core surface can be neglected. The contours 
are depicted with the probability intervals of 0.1, assuming /i^ = 10^^^/is and cos 26* = 1. It can 
be seen that some observable effects on the spectrum of the emitted v^J's are expected if Bq is 
stronger than ~10^ [Gauss] and Am? is larger than ^^10"^ [eV^j. Note that the typical energy 
range of the neutrinos which are observed in a water Cerenkov detector is 10-70 MeV. The lower 
margins of the strong conversion region {P > 0.9, where P is the conversion probability) in the 
contour maps are, in all the four models, contours which runs from the upper left to the lower right 
direction. This is due to the fact that the adiabaticity condition is well satisfied with larger values 
of Am? jEy, as discussed in the previous section. We refer to this marginal region in the contour 
maps as "the continuous deformation region" , because the conversion probability continuously 
decreases with increasing neutrino energy in this region and the spectral deformation is expected 
to be continuous. What is interesting about these maps is that some band-like patterns can be seen 
in the relation of the conversion probability and the value of Aw? jEy. For example, the conversion 
probability in the region of Am^jEy = 5 X 10-4-5 x 10"^ [eVVMeV] in Fig. |ll| is much lower 
than that in the other regions of the map. These patterns come directly from the jumps in the 
matter potential due to the onion-like structure of the isotopic composition in giant stars [ pl| . 
When Aw? I El, is in the above region, the resonance in the model NH4 occurs at the surface of the 
iron core and the interval of Aw? jEy corresponds to the jump of AV at the surface (see Fig. ||). 
Since the matter potential changes suddenly here, very strong magnetic field is necessary for the 
satisfaction of the adiabaticity condition, and consequently the resonant conversion is significantly 
suppressed. We refer to such bands of Aw? jE^ as "the weak adiabaticity band", hereafter. At 
each boundaries of the onion-like structure of massive stars, this weak adiabaticity band appears 
due to the jump in the matter potential. It should also be noted that in the weak adiabaticity 
bands, the conversion probability only weakly depends on Aw? jE^ because the location of the 
resonance layer is not changed in a band. In Figs. |-jl|, we can see the difference between the WW 
and NH models as well as between 15 and 25 M© models. Although there are some quantitative 
differences, almost all qualitative features are the same for these four models. 

Next we show the contour maps of the conversion probability for the zero metallicity models, in 
Figs. |l^ and ^ (for the models WW15Z and WW25Z, respectively), with the region of Am^ /E^ = 
10-"-10-^ [eVVMeV] and Bq = 10^-101°. When Am^/E^ is larger than ~ 10"^ [eVVMeV], 
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the profile of the contour maps is quahtatively similar to that of the solar metallicity models. But 
the adiabaticity condition can be satisfied with smaller strength of magnetic fields and the region 
of complete conversion becomes somewhat larger, because in the inner part of the isotopically 
neutral region (0+Si and 0+Ne+Mg layers) the value of (1 — 2Ye) is about 1 order of magnitude 
smaller than that in the solar metallicity models. Especially, in the model WW25Z, AH and 
ji^B arc comparable in this region (sec Fig. ^ and the precession effect is no longer negligible, 
leading to the more complicated feature of the contour map of WW25Z than of WW15Z. When 
Am'^/Ei, < 10~^ [eV^/MeV], the strong precession occurs in the outer part of the isotopically 
neutral region (0+C and He layers), where fi^B is much higher than AH. In contrast to the 
solar metallicity models, the conversion still occurs with such a low value of Am'^ /E^, even when 
Am^ = 0. Further interesting is that with Aw? /Ei, lower than ~10~^ [eV^/MeV], the conversion 
probability changes periodically with /ii,i?o (Fig- ^3|). This can be understood as follows. The 
precession effect in the outer part of the isotopically neutral region is very profound and then 
this precession is stopped almost suddenly at the boundary of the helium layer and the hydrogen 
envelope where \AV\ increases by 5-10 orders of magnitude. The final phase of the precession 
strongly depends on ^.^Bq, because the precession length is given as 



Therefore the conversion probability changes periodically with ^^Bq. The examples of strong 
precession effect are shown in Fig. ^ using the WW15Z model with Am? = and some values of 
-Bq. The conversion probability as a function of the radius is shown in this figure. The precession 
begins at r = O.O25i?0 and ceases at r = 0.28i?Q (see also Fig. P). One can see that the precession 
length becomes larger with propagation of neutrinos, because B decreases with r. It is also clear 
that the change in Bq leads to the change of the precession length, and hence to the oscillation 
of the final phase of precession. Below Am? / E^ ^ 10^^^ [eV'^/MeV], the effect of the mass term 
can be completely neglected and the conversion probability becomes constant (but never vanishes) 
with neutrino mass or energy. 



B. Spectral Deformation 

All of the qualitative features of the spectral deformation due to the resonant spin-fiavor con- 
version of ^ are clearly understood by the contour maps given in the previous section (Figs. 
^-|lj). The most easily detectable fiavor in a water Cerenkov detector is v^'s because of the large 
cross section of the reaction D^P ne~^ , and they are detectable above the positron energy of ~5 
MeV in the Supcr-Kamiokandc detector, which has the fiducial volume of 22,000 tons If we 
consider the positron energy range of 10-70 MeV, which includes almost all of the events, this 
range corresponds to a vertical bar in the contour maps with fixed values of Am^ and Bq. The 
samples of such a bar are shown in Figs. Ill] and |l2| (NH models), and the corresponding spectral 
deformation of the events at the SK are shown in Fig. ^ The distance of the supernova is set to 
10 kpc and the total energy of each type of neutrinos is assumed to be 5 x 10^^ erg. We use 5 and 
8 MeV as the temperature of j/e's and Z/'^'s, respectively, and the Fermi-Dirac distribution with 
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zero chemical potential is assumed for both D^s and z^^'s p7[ |. The cross section of the dominant 
reaction of PeP — > ne'^ is 9.72 x lO^'^^E'ePe cm^ where and Pe is the energy and momentum 
of recoil positrons. The appropriate detection efficiency curve is also taken into account |^^. The 
thick solid line in Fig. |l^ is the expected differential event number of Ve's without any oscillation 
or conversion. 

As mentioned in the previous section, the three characteristic regions appear in the contour maps 
for the solar metallicity models: A) the complete conversion region, B) the continuous deformation 
region, and C) the weak adiabaticity band. (A, B, and C correspond to those in Figs. |l^ and 
p^ .) When the conversion is complete, we can see the original spectrum as Pg's, and the event 
number is considerably enhanced because of the higher average energy (thin solid line in Fig. p^ . 
When the vertical line in the contour map lies in the continuous deformation region, conversion 
probability decreases with increasing energy of neutrinos and consequently the original v^^s are 
dominant in the lower energy range, while the original De's are dominant in the higher energy 
range (short-dashed line in Fig. |l6|, also corresponding to the vertical line (B) in Fig. 12). Note 
that this feature is based upon the assumption that the radial dependence of the magnetic fields 
is a dipole {B cx r~'^) and fi^B drops faster than \d{AV)/dr\^/'^ with increasing radius. On the 
other hand, in the weak adiabaticity band, the energy dependence of conversion probability is 
rather weak (long-dashed line in Fig. |l6|, also corresponding to the vertical line (C) in Fig. pi] ). 
Because the resonance always occurs in the same place (jumps in the matter potential), this feature 
does not depend on the assumption of the radial dependence of magnetic fields, in contrast to the 
case (B). Finally, quite interesting deformation is expected if the vertical line in the contour maps 
crosses the boundary of the weak adiabaticity band (the vertical line (D) in Fig. 11). Because the 
conversion probability changes almost suddenly at the boundary, the spectrum of the event rate 
suffers drastical deformation at a certain positron energy (dot-dashed fine in Fig. |l^). The used 
values of (Am^ [eV^], Bq [Gauss]) for the vertical lines A, B, C, and D in Figs. |ll| and |l^ are 
(5x10-3,7 X 10^), (5 X 10-4,5 X 10^), (3 x 10-'\5x 10^), and (1.5 x 10^2,2 x 10^), respectively. 

In the zero metallicity models, the feature of the spectrum deformation is similar to that of the 
solar metaUictiy models when Am^/^:^ > 10"'^ [eV^/MeV]. However, if Am^ lO-'^-lO-'' [eV^], 
the conversion probability increases with neutrino energy, because the precession in the outer part 
of the isotopically neutral region becomes effective (Fig. |l^). With Airi^/Ei, < 10-^ [eV^/MeV], 
the conversion probability becomes constant with neutrino energy, but changes periodically with 
Bq- In the model WW25Z, the precession in the inner part of the isotopically neutral region (0+Si 
and 0-|-Ne-|-Mg layers) is also effective, and the probability may change rapidly and complicatedly 
with neutrino energy (Fig [l^ ) . 



V. DISCUSSION 

We found that the difference of the stellar metallicity significantly affects the resonant spin-flavor 
conversion of De ^ i^fi, and some implications from this fact are discussed in the following. The 
lifetime of massive stars which end their life by the gravitational collapses is very shorter than 
that of the Sun, and the progenitors of observed supernovae are therefore younger. Consequently, 
the metallicity of the Galactic supernova is expected to be at least the solar abundance or more 
metal-rich. If the metallicity is higher than that of the Sun, the suppression of (1 — 2Ye) will be 



12 



weaker and the Bq which is required to satisfy the adiabaticity condition becomes larger. On the 
other hand, the Large and Small Magellanic Clouds are known to be very metal-poor systems [ p9[ . 
Therefore, the resonant conversion will occur with smaller magnetic fields in supernovae in the 
Magellanic Clouds, and also the precession effect may be observed. The another object which has 
relation to the metallicity effect is the supernova relic neutrino background (SRN) |40|. Because the 
SRN is the accumulation of neutrinos from supernovae which have ever occurred in the universe, 
the SRN includes neutrinos from supernovae with quite low metallicity in the early phase of galaxy 
formation. The conversion of Pe ^ conpensates the energy degradation due to the cosmological 
redshift effect and enhances the expected event rate of the SRN. However, because of the small 
expected event rate at the SK |4^, it will be difficult to get some decisive information on the 
spin- flavor conversion from the observation at the SK. 

We did not consider the flavor conversion (the MSW effect) in this paper although this can occur 
with appropriate generation mixing and neutrino masses. With the same lS.m? /E^, the resonance 
of the flavor conversion occurs in more outer region than the resonance layer of the spin-flavor 
conversion, because the matter potential for the flavor conversion are not suppressed by (1 — 2Ye) 
[ prf . Even if the flavor conversion occurs in more outer region, the spectrum of is not changed. 
There may be interesting effect if we consider the conversion in the iron core, or the mutual effect 
of spin- flavor and flavor conversion among the three generations of neutrinos, as pointed out by 
the earlier publication 1 1 1 1 . 

The turbulence in the radial dependence of the magnetic fields was ignored in this paper. In the 
solar metallicity models, AiJ is much larger than ji^B except at the resonance, and the conversion 
probability is determined only from the strength of the field at the resonance layer. Therefore 
the turbulence does not affect the evolution of conversion probability of neutrinos. However, the 
turbulence disturbs the relation of Bq and B at the resonance layer. When the neutrino energy 
changes, the location of the resonance layer also changes, and hence the neutrino spectrum can be 
disturbed by some strong turbulence in the magnetic fields. 

The effect of dynamics in the collapse-driven supernova was also not taken into consideration. 
Although it seems unlikely that the shock wave is propagated through the whole isotopically neutral 
region in a few tens of seconds, the inner part of the isotopically neutral region may be dynamically 
disturbed by the shock wave. If Am? /Ei, is large, the resonance occurs at the inner part of the 
isotopically neutral region, and the dynamical effect may change the situation of the resonant 
conversion of the neutrinos emitted in the later phase of emission (> 10 sec after the bounce). We 
give here simple discussion on the effect of the change in density, assuming that Ye is conserved 
during the shock propagation. (The composition of matter is drastically changed by the shock wave, 
but the change in Yg requires the weak interaction.) The matter potential (Ay) changes as oc p. 
On the other hand, if we assume the conservation of the magnetic flux in dynamical plasma, the 
field strength changes as oc p^^'^, and hence the precession becomes more effective with decreasing 
matter density (see Eq. (|^)). However, the adiabaticity condition is the competition of p,^B 
and |d(Ay)/dr|i/2, and if we assume \d{A.V)/dr\ scales as (length) (homogeneous expansion 
or compression), the scaling of \d{/S.V) / dr\^/'^ is the same with that of magnetic fields, cx p^^^. 
Therefore, the adiabaticity condition is not strongly affected by the dynamics of the shock wave. 

It should also be noted that the observed data of neutrinos from SN1987A |28||2g| ] favor a softer 
neutrino spectrum than theoretically plausible spectrum of electron antineutrinos. If j/g's are ex- 
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changed with j/^-hke neutrinos {v^, Vt, and their antiparticles) which have higher average energy, 
this discrepancy becomes larger. From this view point, an upper bound on the conversion prob- 
abihty of Pe's and f^-like neutrinos has been derived: P < 0.35 at the 99 % confidence level 
The earlier paper of the spin-flavor conversion ||l^] used this constraint in order to derive 
an upper bound on the neutrino magnetic moment. However, the above constraint on P and its 
confidence level suffer considerable statistical uncertainty because of the small number of events 
in Kamiokande and 1MB. Therefore, we avoid a decisive conclusion about the upper bound on ii^, 
from SN1987A, although the strong conversion region (P > 0.9) in the contour maps (Figs. |-|l4l) 
may be disfavored. Also we point out here that the metallicity effect on the conversion probability 
should be taken into account when one attempts to constrain fi^, from the SN1987A data, because 
the Large Magellanic Cloud is a low-metal system. 



VI. CONCLUSIONS 

Neutrino spin-flavor conversion of <-!■ i^^ induced by the interaction of a flavor changing 
magnetic moment of Majorana neutrinos and magnetic flelds above the iron core of collapsing 
stars was investigated in detail. The effective matter potential of this conversion mode (Pg t'/i) 
is proportional to (1 — 2Ye), and hence this value is quite important to this resonant conversion. 
However, this value is determined by isotopes which are quite rarely existent, and in order to 
estimate the effect of the astrophysical uncertainties, we used the six precollapse models, changing 
the stellar masses, metallicities, and authors of the models. The components of Hamiltonian in 
the propagation equation (||) are shown for all the models in Figs. and qualitative features of 
the conversion can be understood from these figures. The results of the numerical calculation for 
all the models are shown in Figs. as contour maps of conversion probability as a function of 

the two parameters of /S.rn? /Ei, and i?o, where So is B at the surface of the iron core. 

For the solar metallicity models, observable effects are expected when /S.m? /E^ is in the range of 
IQ-S-lQ-i [eVVMeV] and jji^ > 10-^^{10^G/ Bq) [hb] (Figs. The difference of the stellar 

masses leads to the different thickness and location of the layers of the onion-like structure in 
massive stars and this effect appears in the contour maps, although the effect is rather small. The 
qualitative features of the contour maps for the WW and NH models are also not so different, and 
the model dependence of the conversion probability can be roughly estimated by the comparison 
of these figures. Although the dependence on the stellar models or stellar masses is rather weak as 
shown in Figs. p|-p^, it was found that the metal abundance of the precollapse star significantly 
affect the value of (1 — 2Ye). The difference between the solar and zero metallicity is prominent 
especially in the 0-l-C and He layers, and the strong precession between i^fj. occurs with 

small Am'^/E,^, because ^„B is much larger than Am'^/Ei, in this region. In contrast to the 
solar metallicity models, the conversion occurs even when Am^ = 0. The probability changes 
periodically with Bq because of the precession effect. (See Figs. ^ and [T^ .) 

Considering the above properties, the expected spectral deformation of D^s can be summarized 
as follows. For the solar metallicity models, there are roughly three types of the energy dependence 
of the conversion probability: 1) complete conversion in a range of neutrino energy, 2) conversion 
probability decreases with increasing energy when the energy range is in 'the continuous defor- 
mation region' in the contour maps, and 3) incomplete conversion and weak energy dependence 
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of conversion probability when the energy range is in 'the weak adiabaticity band' in the contour 
maps. Examples of these types of spectral deformation are given in Fig. 16, (For the explanation 
of 'the continuous deformation region' and 'the weak adiabaticity band', see section |^.) Fur- 
thermore, there appear some interesting jumps in the spectrum if the energy range of 10-70 MeV 
includes boundaries of the complete conversion region and the weak adiabaticity band. Irrespec- 
tive of the stellar models or stellar masses, such a boundary exists especially at the surface of the 
iron core, where the matter potential suddenly changes. For the zero metallicity models, although 
the feature of the spectral deformation is similar to that of the solar metallicity models when 
^w? jEy > 10~^ [eV^/MeV], energy-independent conversion is possible with quite small /S.rn? /Ei, 
(< 10^^ [cV^/MeV]), which does not occur in the solar metallicity models. 
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TABLES 



TABLE I. The table of (1 - 2Ye) in the isotopically neutral region (O+Si, O+Ne+Mg, 0+C, and He 
layers) for all the models of precollapse stars used in this paper. The models WW15 and WW25 are 15 
and 25 Mq models by Woosley and Weaver (1995) |25|. There are models with two different metallicities 
in the WW models, and they are distinguished by 'S' (solar metallicity) and 'Z' (zero metallicity) . The 
models NH4 and NHS are 4 and 8 Mq helium core models (corresponding approximately to 15 and 25 
Mq main sequence stellar masses) by Nomoto and Hashimoto (1988) [Q. The isotopes in the parentheses 
under the values of (1 — 2Ye) are the main isotopes which cause the deviation of Ye from 0.5. 

Model 0+Si O+Ne+Mg 0+C He 

WW15S 2 xlO"^ 2 xlO"^ 2 xlO"^ 1-10 xlO"* 

(3*S, ^«Ar) (^^Na, 25,26^g) (22Ne, ^^Mg) (^^O, ^^Fe) 

WW25S 2 xlO"^ 2 xlO"^ 2 xlO"^ 2 xlO"^ 

(^%^»Ar) (^^'^^Mg^^^Al) (^^Ne, 25'26Mg) (^^Ne) 

WW15Z 6-10 xlO"" 4 xlO"* 1-6 xlO"* 2-100 xlO"^° 

(27^1^ 29,30gi) (23Na, 26Mg) (^^Ne, ^^Mg) (i80) 

WW25Z 4 xlO"'' 1 xlO"* 2 xlO"'' 2 xlO"'' 

(2^A1, »"Si) (^^Na, 26j^g) (^^Ne, 25'26Mg) (i^.^C, ^^Ne) 

NH4 1 xlO"^' 2 xlO"^ 2 xlO"^ 2 xlO"^ 

C% ^»Ar) (25,26^g^ 27^1) (22j^g^ 26^g) (18q^ 22j^g) 

NHS 2 xlO-^' 1 xlO"^ 1 xlO"^ 2 xlO"^ 

(27^1_ 29,30si) (23Na, 2^-26Mg) (^^Ne, ^^Mg) (^^O, ^^Ne) 
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Figure Captions 



FIG. 1. The components of the Hamihonian in the propagation equation (g) are shown for the Woosley 
& Weaver's solar-metalhcity 15 Mq model, as a function of the radius from the center of the star. The 
absolute value of the effective matter potential, AV, is shown by the thick, solid line. The absolute value 
of the difference of the diagonal components of the Hamiltonian, AH, is shown by the dashed lines for 
the two values of the mass term, Arvi^ jEy = 10^'' and 10^^ [eV^/MeV]. The resonant conversion occurs 
where AH — (the resonance layer). The dotted lines show fj,^B, i.e., the off-diagonal components of the 
Hamiltonian. The indicated values of B are those at the surface of the iron core (Bo) in units of Gauss, 
assuming /i^ = 10^^^/is. The thin solid line shows \d{AV) /dr\^'''^ [eV], and the adiabaticity condition 
is satisfied when ^^B > \d{AV) /dr\^'''^ at the resonance layer. The dot-dashed line shows ^^Bth, again 
assuming = 10 fiB, where Bth is the magnetic field whose energy density is the same with the thermal 
energy density in the star. The dominantly existent nuclei in each layer of the precollapse model are shown 
in the top of this figure. 

FIG. 2. The same as Fig. |l|, but for the Woosley & Weaver's solar- met allicity 25 Mq model. 

FIG. 3. The same as Fig. |l|, but for the Nomoto & Hashimoto's 4 Mq helium-core model. 

FIG. 4. The same as Fig. |l|, but for the Nomoto & Hashimoto's 8 Mq helium-core model. The 0-|-Si 
and 0-fC layers are negligibly thin in this figure. 

FIG. 5. The same as Fig. |l|, but for the Woosley & Weaver's zero-metallicity 15 Mq model. Note that 
the scale of the vertical axis is different from the previous four figures. 

FIG. 6. The same as Fig. 0, but for the Woosley & Weaver's zero-metallicity 25 Mq model. 

FIG. 7. The evolution of conversion probability with the radius from the center of the star. The 
Nomoto & Hashimoto's 4 Mq helium-core model is used, and Air? j Ey is assumed to be 10~''[eV'^/MeV]. 
The four different values of So are used for the calculation, where Bo is the magnetic field at the surface 
of the iron core. The neutrino magnetic moment, p^v, is assumed to be 10~^^/is. 

FIG. 8. The same as Fig. |, but for the case with Arri^ / = 10"^[eV^/MeV]. The three values of Bo 
are used as indicated in the figure. 

FIG. 9. The contour map of conversion probability for the Woosley & Weaver's solar-metallicity 15 Mq 
model, as a function of Am^ jEy and Bo, where Bo is the magnetic field at the surface of the iron core. The 
neutrino magnetic moment, /i^, is assumed to be XQT^^ [ib- The contours are depicted with the interval of 
0.1. 
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FIG. 10. The same as Fig. H, but for the Woosley & Weaver's solar-metaUicity 25 Mq modeL 



FIG. 11. The same as Fig. but for the Nomoto & Hashimoto's 4 Mq helium-core model. The two 
vertical bars, (C) and (D) correspond to the energy range of 10-70 MeV with (Am^ [eV% Bo [Gauss]) = 
(3 X 10~^, 5 X 10®) and (1.5 x 10"'^, 2 x 10^), respectively, where the energy range is that of positrons 
which are observed in a water Cerenkov detector. 

FIG. 12. The same as Fig. ^ but for the Nomoto & Hashimoto's 8 Mq helium-core model. The values 
of (Am^ [eV^], Bo [Gauss]) for the two vertical bars, (A) and (B) are (5 x 10"^, 7 x and (5 x 10"*, 
5 X 10®), respectively. 

FIG. 13. The same as Fig. ^, but for the Woosley & Weaver's zero-metallicity 15 Mq model. Note that 
the region of Am^/i?^ is different from that in Figs. p|-|l^ 



FIG. 14. The same as Figure |l^, but for the Woosley & Weaver's zero-metallicity 25 Mq model. 

FIG. 15. The same as Fig. ^, but for the Woosley & Weaver's zero-metallicity 15 Mq model and 
^^rr? jEv = 0. The three values of Bo are used as indicated in the figure. Strong precession occurs in the 
region of 0.025-0.28 Rq. 

FIG. 16. Examples of the spectral deformation of j/^^s for the solar metallicity models. The expected 
time-integrated differential events at the SK detector are shown, assuming that the supernova occurs in the 
Galactic center {D — 10 kpc). The upper box shows the corresponding conversion probability of Ve ^ ffi- 
The thick solid line shows the ordinary expected events of />e's without any oscillation or conversion, and 
the thin solid line the differential events when all of the original i^^'s are completely converted to v^s in 
the whole energy range. The short dashed line is the events when (Am^[eV^], Bo[Gauss]) — (5 xlO~*, 
5 X 10®), using the model NH8, corresponding to the case (B) in Fig. ^ The long-dashed and dot-dashed 
lines are differential events with (Am^, Bo) = (3 x 10~^, 5 x 10®) and (1.5 x 10"^, 2 x 10®), respectively, 
and the model NH4 is used for both the lines. The former corresponds to the bar (C) and the latter to 
the bar (D) in Fig. 0. 
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